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Abstract 

The use of bremsstrahlung photons produced by a linac to induce photonuclear reactions is wide 
spread. However, using a clinical linac to produce the photons is a new concept. We aimed to 
induce photonuclear reactions on zinc isotopes and measure the subsequent transition energies and 
half-lives. For this purpose, a bremsstrahlung photon beam of 18 MeV endpoint energy produced 
by the Philips SLI-25 linac has been used. The subsequent decay has been measured with a well- 
shielded single HPGe detector. The results obtained for transition energies are in good agreement 
with the literature data and in many cases surpass these in accuracy. For the half-lives, we are in 
agreement with the literature data, but do not achieve their precision. The obtained accuracy for 
the transition energies show what is achievable in an experiment such as ours. We demonstrate the 
usefulness and benefits of employing clinical linacs for nuclear physics experiments. 

PACS numbers: 25.20.-x, 23.20.Lv, 27.50.+e, 29.20.Ej 


* haris@akdeniz.edu.tr; http://nukleer.akdeniz.edu.tr/en 
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I. INTRODUCTION 


The study of photon-induced nuclear reactions has been of interest over the years. The 
motivation for studying the interaction of photons with the nucleus are many, ranging from 
the fundamental nuclear structure studies to studying the inner working of a nuclear reac¬ 
tor and the inner processes in a star. Concerning nuclear reactors, photonuclear reactions 
contribute to the general performance of the reactor and their accurate understanding is a 
necessary part of any reactor core simulation. At the same time, the issues of nuclear astro¬ 
physics involve many reactions relating to photons, since, in a star, they are ubiquitous. In 
fact, photonuclear reactions are crucial steps in many of the nucleosynthesis processes gener¬ 
ating the observed abundances of elements. Thus, many photonuclear reaction experiments 
are aimed at and motivated by astrophysical references to nucleosynthesis, see Ref. 
and references there-in. 

In addition, the photonuclear reactions have been frequently applied in photo-activation 
studies pertaining to various fields. This combination of basic research aimed at fundamental 
physics and applied research makes the photonuclear reaction quite interesting and facilitates 
active research in several fields. This activity is also on the increase as the availability of 
radiation sources increases and their cost decreases. For the study presented in this paper, 
the availability of a good radiation source was the main impetus for embarking on this 
research activity. 


Photo-activation experiments involvi^ linear accelerators have been performed at sev- 

inost notably at specialized laboratories such 

iflLi 


eral institutions around the world Ref. M 


as S-DALINAC at 
in Rossendorf Ref. 


"U Darmstadt Ref. llfl lll| and ELBE in Forschungszentrum Dresden 
12| . However, a study in Ref. 13| has shown how a non-specialized 


instrumentation such as a clinical linac can also be employed for such experiments. This 
idea is a departure from conventional approaches in which linear accelerators used for pho¬ 
tonuclear experiments have been designed and commissioned with sole use in nuclear physics 
experiments in mind. 

The concept of using an ’’off-the-shelf” linac, originally designed for a different purpose, is 
a novel idea. Although this idea is an attractive one, the original study presented in Ref. [13| 
has not been followed up on and its impact has been very limited. In light of availability 
of clinical linacs - in Turkey alone there are over 200 such devices - it is clearly a wasted 
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opportunity not to use them for research. Especially important is the potential contribution 
such devices can have on progress of experimental nuclear physics in developing countries. In 
developing countries access to specihcally designed linacs is practically non-existent, whereas 
access to clinical linacs is often readily available. Thus either trying to obtain one such linac, 
after its duty cycle in medicine is over, or working closely with the local hospital, to utilize 
the linac still in us in its off-hours, presents a way how these device can help developing 
countries to contribute to the global nuclear physics knowledge development. It is important 
to note further that there has been a great expansion of available radiation sources in the 
form of clinical linacs not just in developing countries, but across the world. 


The main aim of this work is to follow up on the original work of Ref. 


13| and revive 


this line of research. We will show how a repurposed clinical linac can produce solid nuclear 
physics experiments, and that it is possible even without a sophisticated specifically designed 
linac to obtain good quality data. What we intend to demonstrate is that nuclear physics 
experiments with a clinical linac are a viable concept and one that can be of value to the 
global science community. 

For this purpose it is essential to choose an appropriate example where the opportunity for 
contribution can be clearly demonstrated. A kind of niche where limited budget experiments 
can contribute to the global knowledge. One such opportunity presents itself in cases where 
the measurements of basic nuclear quantities, such as transition energies and half-lives, have 
not been performed for several decades even in the case of elements near to the valley of 
stability. Hopefully by presenting an appealing and convincing example the viability of the 
concept will be clearly demonstrated contributing to its popularity. 

In order to best illustrate the potential of clinical linacs in experimental nuclear physics, 
we have chosen to focus on zinc isotopes to study the transition energies and half-lives of 
isotopes created by photonuclear reactions. The choice of zinc as a target of bremsstrahlung 
photons produced by the clinical linac is motivated by several reasons. Mainly, we are 
interested in studying intermediate mass nuclei, especially those isotopes on the proton-rich 
side of the nuclear chart. In addition, the data for the ®^Zn isotope, the main focus of 


this study, and its /3-decay product ®^Cu are nearly 40 years old, see Ref. 


Q- 


This is the 


second part of the motivation for this study, i.e. to revisit some of the experimental results 
obtained quite some time ago and not investigated again. The goal behind this choice is to 
improve the data accuracy thus illustrating what this and similar studies can accomplish. 
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This work aims at demonstrating both the feasibility of photonuclear experiments with a 


clinical linac as well as to point out that there are isotopes, close to the val 
whose reexamination is worthwhile. In the initial study presented in Ref. 


ey of stability, 


13| the second 


part of this motivation was not explored in detail. In addition, we would like to note that we 
have already had reasonable success with this approach in our preliminary study published 
in Ref. 151. 


The paper is organized as follows; in Section [TTl we give the properties of the clinical linac 
used in this study. Detector setup and the experimental procedure used are presented in 
Section m The explanation of the data analysis and the presentation of the results are 
provided in Section [I^ and hnally, we present our summary and conclusions in Section |Vl 


II. PROPERTIES OF THE CLINICAL LINAC 



FIG. 1. Philips SLI-25 clinical linear electron accelerator of Elekta TM Synergy TM. 


The photonuclear reactions by their very nature require a source of photons, energetic 
enough to excite the target nuclei. In our experiment, the photon source was a clinical 
electron linac. The linac was re-commissioned from its usual role in medical physics for the 
use in the study of photonuclear reactions. It should be noted that a pioneering study by 
Ref. demonstrated that these clinical linacs (cLINAC) have properties appropriate for 
use in photonuclear reactions. In fact, it was demonstrated that they even compare well 
with the linacs made specihcally for physics experiments. And while our study has been 
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conducted with a re-commissioned linac the study in Ref. was done while the linac was 
still in medical use. While the later arrangement is possible with some coordination and 
planing the former arrangement gives more freedom for conduction experiments. 

In our setup, we have used a cLINAC SLI-25 manufactured by Philips Medical Systems 
(currently part of Elekta TM Synergy TM), as a bremsstrahlung photon source, shown in 
Fig. [U The accelerator’s technical documentation can be found in Ref. Ig. The cLINAC 
primary electron beam is generated by an electron gun with an energy of about 50 keV. The 
electron gun in SLI-25 is a diode design with a 400 Hz pulse repetition frequency. After 
injection into the linac’s copper cavity, the electrons are accelerated by a radio-frequency 
wave with 3 GHz (2856 MHz), S-band. The copper cavity is a traveling wave design, where 
the power is injected at the beginning of the accelerator structure. Even though the linac is 
designed to operate at energies up to 25 MeV, the power is provided by the magnetron instead 
of the klystron, which is more common for such energies. The nominal power supplied by 
the magnetron is 2.5 MW at 4 MeV (low energies) and 5 MW at 25 MeV (high energies). 
As is usual, all steering and focusing of the beam is achieved by standard magnetic and 
electrostatic devices. 


Target W (d=0.3 mm) 



sample position-- - ' y 

without cLINAC modification 


electron target 

primary collimator 

flattening filler 
ionization chamber 


Y diaphragm (with MLC) 

I X diaphragm 
■ - Sample 56 cm from target 


FIG. 2. Schematic view of the cLINAC. 


After exciting the copper cavity, the beam is steered through 112.5° slalom-type magnetic 
and electrostatic devices before falling onto a high-Z element target, in our case, tungsten. 
The tungsten target is 0.3 mm thick and serves as an electron stopper and a bremsstrahlung 
photon source. Because the linac’s original purpose was medical, i.e. for patient treatment. 
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the bremsstrahlung photons are collimated and flattened with several Alters placed as shown 
in Fig. [21 The resulting photon beam is spatially uniform with no position dependence, i.e. 
the beam has the same intensity at the center as it has at the edges of the field. To further 
illustrate this, we show in Fig. [3] a simulation of the photon flux created by the linac. As 
is evident in the figure the photon flux is flat across the sample position at 56 cm. For 
us, this is especially convenient as it avoids the necessity for complicated arrangements to 
create uniform sample irradiation. The focusing and collimation is a standard feature of 
all cLINACs, as it is paramount to maintain excellent spatial dose profile control when 
irradiating a patient. In fact, it is common to require dose knowledge to better than 3% 
accuracy. The dose and the dose spatial distribution are measured regularly as part of the 
standard linac performance monitoring. 


Linac simulation (Normalized flux) 



z [cm] 


FIG. 3. Linac photon-flux distribution simulation normalized to the maximal flux. 


Fig. [3] shows a simulation of the photon flux produced by the linac. The figure shows the 
photon flux relative to the maximum flux in arbitrary units. It shows the location of the 
collimators and the flattening Alter as well as their influence on the flux. Once the beam 
passes the flattening Alter, which is of conic shape and much thicker in the center than at the 
edges, the beam loses any spatial dependence in intensity and becomes flat across the width 
of the Alter. In addition, the collimators strongly reduce the beam intensity outside the focal 
area leading to a difference of several orders of magnitude in intensity inside and outside of 
the focal opening. After the flattening Alter and the collimators the beam remains both flat 
and focused as it continues its trajectory. Once it reaches the sample target position at 56 
cm, an area of approximately 20x20 cm^, it still has a uniform beam profile. The scoring 
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plane for the bremsstrahlung spectrum was placed at the surface of the water phantom at 
100 cm downstream. At this distance the beam covers uniformly an area of 40x40 cm^. 
Note that the surface of the water phantom was greater than 40x40 cm^. Even as it travels 
through the phantom the beam still retains its uniformity to a very good degree. Although 
not used in the present experiments the collimators and the flattening hlters can be used to 
create narrower or irregularly-shaped held shapes. 


Under all of these conditions, the question of bremsstrahlung energy distribution, which 
comes from the linac, becomes a quite complicated one due to presence of all the hlters 
and collimators. A simulation of the photon energy distribution coming from the SLI-25 


was performed using the BEAMnrc package. Ref. 


17| . The resulting distribution is shown 


in Fig. m The distribution shown was calculated with an electron beam accelerated over 
an 18 MV potential diherence impacting on a tungsten target of 0.3 mm thickness. As 
is usual, the distribution was calculated per incoming electron. For the estimation of the 
full hux, it is necessary to take into account the number of incident electrons. In our case, 
this number can be estimated from the dose delivered by the beam, 5 Gy/min Ref. Ig, to 
be about 10^^ electron/s. Combined with the simulation from BEAMnrc, this gives about 
5 X 10^photons/(MeVcm^s) at E = 6 MeV. 



FIG. 4. EGS simulation of bremsstrahlung energy distribution from an electron beam accelerated 
over an 18 MV potential difference impacting on thick tungsten target in Philips SLI-25 clinical 
linac. The distribution shown is per incident electron. 





III. DETECTOR SET-UP AND EXPERIMENTAL PROCEDURE 


In the experiment, one zinc sample was placed immediately outside the cLINAC head at 
about 56 cm source to sample distance, see Fig. [2l The sample in question was a zinc disk 
weighing 5 g with 5 mm diameter and 1 mm thickness. The sample was irradiated for about 
35 min. The time was optimized such that a good count rate was achieved in the detector; 
the count rate was low enough not to cause pile up in the detector and yet high enough to 
observe several relatively weak peaks (intensity ~ 0.05%). The sample was made primarily 
of zinc (67% by mass fraction), but it had several impurities, primarily magnesium, some 
aluminum, little iron and trace amounts of other elements. However, in the previous tests, 
we have already observed that the impurities did not show up in the spectrum. In addition, 
since no special effort was made to isotopically enrich the target, the presence of different 
zinc isotopes can be considered as those of their natural abundances (®^Zn - 49.17%, ®®Zn - 
27.17%, ®^Zn - 4.04%, ^^Zn - 18.45% and ^°Zn - 0.61%). 


After the irradiation, the sample was transported to the Physics Department of Akdeniz 
University, where the detector was located. The detector used is a high-purity Germanium 
detector (HPGe). It is a p-type, coaxial, electrically cooled HPGe detector, placed in a 
well-shielded cavity. The shield is 10 cm thick lead with an inner surface covered by a 2 mm 
copper foil to reduce the Pb X-rays generated in the Pb shielding. The HPGe detector used 
is a gamma-ray spectrometer from AMETEK-ORTEG (GEM40P4-83) with 40% relative 
efficiency and resolution of 768 eV FWHM at 122 keV for ^^Go source and 1.85 keV FWHM 
at 1332 keV for ®°Go source Ref. Q . It is connected to a set of Nuclear Instrumentation 
Modules consisting of ORTEG preamplifier, bias supply, spectroscopy amplifier, analog-to- 
digital converter and a computer. Data acquisition was carried out with MAESTR032 
software. Through some trial and error in preliminary studies with zinc, we came to the 
conclusion that the best setup of electronics and detector for this experiment is the one with 
16830 channels and a channel step of about 0.18 keV/channel, giving about 3000 keV for 
the last channel and the maximum transition energy which can be observed. 


The sample was placed in front of the detector about 10 min after the irradiation and 
7 -ray counting continued for three days. During the counting, spectra were automatically 
recorded at regular time intervals. Initially, those time intervals were short ~3 s, designed 
to follow the short-lived isotopes, while the later ones became longer ~20 min when focusing 
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on longer-lived isotopes. 

Immediately before counting, a set of calibration sources were measured. For calibration, 
we used two sets of sources. The point sources were supplied by the Qekmece Nuclear Re¬ 
search and Training Center (IAEA 1364-43-2) and contained Co-60, Na-22, Mn-54, Cd-109, 
Co-57, Cs-137 and Ba-133 isotopes. The second set is a soil sample provided by the Turk¬ 
ish Atomic Energy Authority (TAEK). It contained different natural radioactive isotopes 
(40K,^26j^g^ and 232xh) with known activities. From this sample, we used only the strongest 
peaks for calibration. 

After the 7 -ray measurement of the sample was performed, an equivalently long natural 
background spectrum was recorded. Once the background spectrum was recorded, the 
experiment was concluded with a second measurement of the calibration sources. The aim 
of this second calibration-sources measurement was to check the stability of the electronics 
and the measurements. In this way, we were able to track any channel shift during the 
measurement and also, through combining the before and after calibration results, eliminate 
or at the very least reduce any systematic errors coming from the channel shift. 

The process of photo-activation applies to all the nuclei for which the bremsstrahlung 
radiation exceeds the reaction threshold. In a way, it is reasonable to assume that all the 
processes that are physically possible are realized inside the sample. The neutron separation 
energies of the zinc isotopes are in the range from 7 to 12 MeV and the proton separation 
energies are in the range from 7.5 to 11.5 MeV. Given that the endpoint energy of our beam 
was 18 MeV, it is expected that all of the stable zinc isotopes will get activated. However, due 
to limitations of the experimental setup, not all 7 -decays from zinc isotopes were observable 
or have been observed. The photonuclear reactions whose signatures we have observed in 
the spectrum are: 


®^Zn + 7 Zn + n. 

( 1 ) 

®®Zn + 7 ^®'^Zn + n. 

( 2 ) 

™Zn + 7 Zn* + n, 

( 3 ) 

®®Zn + 7 Cu + p. 

( 4 ) 


The most important limitation for this reduced list is the half-life of the created nuclei. 
It is because of this that we do not see evidence of ®^Zn since the half-life of ®®Cu is only 5.12 
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min. At the same time the half-life criteria works in our favor and we are able to observe a 
long-lived (13.76 h) isomeric state of ®®Zn. 

Our experimental setup was not designed for direct observation of the above reactions, 
i.e. prompt gammas, thus we have to rely on the radioactive decay of the created nuclei 
and their end-stage 7 -decays. The experiment was intended to observe the half-life of the 
radioactive nuclei produced and the transition energies of their daughter products, not to 
observe the levels of stable nuclei which are used as the targets. In this respect, we list the 
decay reactions which we studied 


®^Zn 

Gu* + e+ + z/. 

( 5 ) 

6 ®Zn 

Gu* + e+ + z/. 

( 6 ) 

69zn* 

^^''Zn + 7 , 

( 7 ) 

®^Gu 

Zn* + e + i). 

( 8 ) 


Here too, we have yet another limitation as to what we can observe. As is noticeable the 
/3-decay of ®®Zn is missing from the above list since it was not directly present in the ob¬ 
served gamma spectrum. The observational limitation here is the low branching ratio of the 
®®Zn /1-decay into the excited states of ®®Ga. Namely this /1-decay goes, almost exclusively 
(99.998%), to the ground state of ®®Ga leaving no gamma transition to be observed by the 
detector. 

In addition to the process of photonuclear reactions, explained above, it should be noted 
that there exists a small ambiguity in the reactions for a few of the cases presented. This 
ambiguity appears due to somewhat high endpoint energy of the bremsstrahlung spectrum. 
As can be noticed, at 18 MeV the photon energy is higher than the neutron separation en¬ 
ergy for most elements. As such, it is expected that the radiation source and the collimation 
materials will become sources of neutrons. Therefore, some of the nuclei observed could 
have been produced in neutron-capture reactions. However, this source of neutrons as sec¬ 
ondary particles is rather low and the processes of photonuclear reactions presented above 
are dominant. Therefore, for the results on which this work focuses, such small interferences 
are not of significance. Only for the cross-section determination can such issues play a role. 

These neutrons coming from the radiator and the collimation materials in the head of the 
linac are of interest in medical therapy and can be measured. In our case we have measured 
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the slow neutron flux as being close to 8000 neutrons/(cm^s) at 100 cm distance from the 
radiator. The measurement was performed in the center of the photon beam. However it is 
expected that the neutrons will have a similar presence even outside of the photon beam, 
since the radiator emits neutrons in all directions. In medical therapy with photons this 
presents a problem since it causes a non-negligible dose outside the intended area. Other 
treatment modalities with charged particle offer better dose profiles and less neutron leakage 


Ref. 


19|, 


2Q|. 


IV. DATA ANALYSIS AND RESULTS 


The data acquisition was performed with MAESTRO software; however, due to its lim¬ 
itations peak analysis and energy calibration were performed with other programs. Peak 
analysis was performed with RadWare code developed for analysis of gamma-ray coincidence 
data, by David Radford of the Physics Division at Oak Ridge National Laboratory Ref. 


while energy calibration was performed in ROOT developed by CERN group Ref. 


m, 


22|. The 


motivation for using two different programs was to strike a balance between the desired 
accuracy and time consumption for the data analysis. At the end, the data were combined 
in a simple calculation sheet, giving the transition-energy value and the associated error. 

Thus, all of the recorded sample and background spectra as well as the calibration spectra 
taken before and after the experiment were analyzed with the RadWare package. The 
strength of the RadWare package and its suitability for the analysis of spectra recorded 
with a HPGe detector lies in the fact that it fits a Gaussian, a skewed Gaussian, and a 
smoothed step function to any number of chosen peaks. In this way, we determine the 
centroids, areas and respective statistical errors of all peaks found in a spectrum. 

The first step of the analysis was performed on the calibration data taken before and 
after the experiment. As is usual, the peak position in terms of channel, i.e. centroid, was 
obtained. The two calibrations were then combined taking into account both the errors 
obtained during the fit as well as the distance between them. Ref. 23j. In this way, we 


obtained a unified determination of calibration-source peak centroids. These were paired 
with the corresponding transition-energy data taken from literature. The best example of 
this procedure is the Go-60 with its two peaks at 1173.237±0.004 keV and 1332.501±0.005 
keV. The centroids and energies are fitted with the aid of ROOT and the energy calibra- 


12 








tion with accurate uncertainty for the fit parameters is obtained. In addition, the fitting 
subroutine provides the correlation between the fit parameters. We note that the ROOT 
fitting procedure takes into account both the uncertainties in energy as well as the centroid 
according to the effective variance method. 

At this point, the question arises as to which is the appropriate fitting function for the 
energy calibration. Nominally, the energy calibration should be linear, but in practice it 
may not be. Thus, we tested linear, quadratic and cubic energy calibrations. In our case, 
the quality of the fit improved as we used higher order polynomials. In fact, the y^/n.d.f. 
for linear fit was 15.48, for quadratic 5.65 and cubic 1.41. Hence, we used the cubic fit for 
the energy calibration. 

The propagation of error from the defining energy calibration equation E = 
then takes all errors into account. Error of the fit parameters cXa., covariance covij or corre¬ 
lation coTij matrix and the errors of the centroid determination itself ach- The error formula 
can then be written as: 



With the calibration done, we proceed to determine the peak position in the measured 
zinc sample spectrum. In addition to the sample spectrum, we also determined the recorded 
peak position in the background spectrum and used it to discard the same peaks found in the 
sample spectrum. What remained was the centroid position information for the peaks, which 
could be assigned to the sample. The energy and uncertainty were calculated according to 
Eq. ffTOj) and are listed in Table [H 

When we compared the energy results, from Table [H with the literature values, it was 
possible to assign these peaks to q-ray transitions in specific isotopes. In Fig. |51 we show the 
irradiated-sample spectrum without any background subtraction. The energy assignment of 
peaks is shown in the figure. All the unassigned peaks are either background peaks or sum 
and escape peaks. On the other hand, all the assigned peaks can be connected with the 
nuclei produced in the photonuclear reactions on zinc isotopes. Hence, it is safe to say that 
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any elemental pollution we might have had in the sample did not show up in the spectrum 
in the form of an unidentihable peak. In this way, we verify that our assumption of de facto 
working with a pure zinc sample was justihed. 
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FIG. 5. Zinc sample spectrum with no subtractions, for the irradiated sample after 3 days of 
counting. Assigned peaks are labeled with their energy in keV based on the energy calibration, and 
also with the decaying isotope. 


Finally, the results obtained for the average energy and the combined variance are com¬ 
pared to literature values in Table [B The literature results quoted are taken from NUDAT 


and come from Nuclear data sheets publications, which for these elements are Ref . 
In turn Nuclear da 


decay. Ref. 


28 


a sheets results are taken from individual publications: Ref. 


3-27|. 


14| for ®^Zn 


29 1 for ®^Cu decay. Ref. 30, M| for ®®Zn decay, and Ref. 32| and others 


for ®®Zn decay. In order to illustrate the quality of the agreement, we also show in Table H] 
how distant are our values from the literature ones. In addition, we show the ratio of the 
distance and the combination of our and literature uncertainties. It may be observed that 
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Ele. 

EAr(keV) 

CJat 

E(keV) 

O'E 

A— \E]\f-E\ 

+ ^e) 

A/D 

63 Zn 

449.93 

0.05 

449.942 

0.017 

0.012 

0.05 

0.24 

63 Zn 

669.62 

0.05 

669.69 

0.03 

0.07 

0.06 

1.17 

63 Zn 

742.25 

0.10 

742.18 

0.03 

0.07 

0.10 

0.70 

63 Zn 

962.06 

0.04 

962.09 

0.04 

0.03 

0.06 

0.50 

63 Zn 

1123.72 

0.07 

1123.80 

0.07 

0.08 

0.10 

0.80 

63 Zn 

1327.03 

0.08 

1327.11 

0.10 

0.08 

0.13 

0.62 

63 Zn 

1374.47 

0.13 

1374.41 

0.08 

0.06 

0.15 

0.40 

63 Zn 

1392.55 

0.08 

1392.56 

0.08 

0.01 

0.11 

0.09 

63 Zn 

1412.08 

0.05 

1412.16 

0.11 

0.08 

0.12 

0.67 

63 Zn 

1547.04 

0.06 

1547.16 

0.12 

0.12 

0.13 

0.92 

63 Zn 

2026.8 

0.3 

2026.70 

0.18 

0.10 

0.4 

0.25 

63 Zn 

2336.5 

0.3 

2336.4 

0.2 

0.1 

0.4 

0.25 

63 Zn 

2536.0 

0.3 

2535.9 

0.3 

0.1 

0.4 

0.25 

63 Zn 

2696.6 

0.3 

2696.5 

0.4 

0.1 

0.5 

0.20 

67Cu 

184.577 

0.010 

184.63 

0.06 

0.05 

0.06 

0.83 

69 Zn 

438.634 

0.018 

438.601 

0.017 

0.033 

0.025 

1.32 

66 Zn 

1115.539 

0.002 

1115.51 

0.06 

0.03 

0.06 

0.50 


TABLE 1. Gamma-ray energies obtained in the present measurement by averaging the results of 
before and after calibration measurements compared to values found in the literature (NUDAT). 

the agreement is good. All the results agree within 1.5a and most are within a or better. 
Even if the uncertainties are not combined, but simply the larger of the two is taken, the 
agreement is still better than 2a. At the same time, the results presented are of similar 
quality or better than the literature ones. The only exception being ®^Zn which is often 
used as a calibration standard and is measured to a very high accuracy. 

In addition to the transition energies, the analysis performed provided information about 
counts and their time evolution. This information can be used to determine the half-life of 
the parent nucleus as the decay of the daughter levels are in secular equilibrium with the 
/3-decay of the parent nucleus. Secular equilibrium is a very good assumption here since the 
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half-lives of the states in question are orders of magnitude smaller than those of the parents. 
The only exception is the ®®Zn, which is an isomeric transition and is thus observed directly. 

Usually the measurement of the half-life has involved the measurement of the decay 
as function of time and the £t to the activity with the exponential decay curve {A{t) = 
AQexp{—Xt)). Also, a naive way would involve hts to the integral of the activity which 
is directly represented by the counts. However, htting the peak in each consecutive step 
independently and then obtaining the half-life from these data is not appropriate since all 
the errors in each step are correlated, but the correlation is unknown. A better approach is 
to integrate the activity in equal-size time steps: 

^T+AT 

C'(T) = / A{t)dt = Coe-^^ (11) 

Jt-AT 

where Cq = Aq/X and T is counting time. So long as the AT is a constant the function only 
depends on T exponentially just like activity. 

In practical terms, this can be performed by taking independent spectra of length AT 
and restarting the count at the end of each time step. In this way, the counts obtained in 
two consecutive spectra are no longer correlated and the errors of the counts obtained in 
this way are independent of each other. Care should be taken to correct the counts obtained 
for dead time since while the detector is busy the sample is still decaying. In our case, this 
was a minor correction since only at the very beginning we had a small amount (~ 2%) of 
dead time. 

In order to simplify the £t, logarithm of Eq. ffTTl) was used for htting. In this case, one 
has a simple linear ht from which the decay constant A is obtained. Once the value of the 
decay constant A has been obtained, it is a straightforward matter to calculate the half-life 
Ti /2 = In 2/A. 

Several examples of the counts dependence on time are shown in Fig. |6l Only the best 
are shown, i.e. the strongest peaks, as well as the dependence of the sum count of all peaks. 
As could have been expected, a linear trend is quite evident and the error-bars are getting 
smaller for stronger peaks and earlier counting intervals. The remaining peaks were too weak 
to be htted, i.e. the htting procedure did not converge or gave very large uncertainties. 

The summary of the hts to the function dehned by Eq. ffTTl) and shown in Fig. [6] is given 
in Table [TTl All of the errors listed in column 3 of Table mi are statistical only. As can 
be seen, the combined-counts ht ohers the smallest error, but the two strongest peaks 669 
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FIG. 6. 


Logarithmic time evolution of counts, as defined by Eq. dm), for all transitions 


combined and separately for a few of the strongest ones. 


Transition 

Ti /2 [min] 

^ stat. 

X^/n.d.f. 

Combined 

38.84 

0.15 

0.782 

449 

42.1 

3.3 

0.700 

669 

38.70 

0.19 

1.616 

962 

39.19 

0.26 

1.003 

1412 

39.2 

1.2 

0.977 


TABLE 11. ®^Zn determined half-life from several transitions. These should be compared with 
NUDAT value of 38.47 ± 0.05 min. 

and 962 are also quite good. When compared to the literature value of 38.47 ± 0.05 min, 
the combined-counts fit is quite decent. The remaining fits scatter noticeably around these 
values. However, all our values are systematically higher than the literature value. This 
systematically higher result is likely the product of isotopic impurity of the sample, i.e. the 
counts of a single isotope are distorted by the presence of other ones. Unfortunately in the 
current experimental setup we can not eliminate this systematic uncertainty. Nevertheless, 
even with these pollution effects, the results are reasonable enough. 

In addition to ®^Zn, we also fitted the y-decay of ®®Zn isomeric transition and obtained 
13.76±0.18 h while the literature value is 13.76±0.02 h in surprisingly good agreement with 
each other. Of the other two observed /3-decays ®®Zn and ®^Cu, we were not able to fit their 
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decay curves. The ®®Zn was not fit due to unfavorable ratio of its long half-life (244 d) and 
observation time (3 d). The ®^Cu was not £t due to the weakness of the signal. Like in 
the case of ®^Zn peaks, not listed in Table [TTl h had too big uncertainties to allow for a 
reasonable £t. 


V. SUMMARY AND CONCLUSION 


In the experiment presented in this paper, we have investigated the spectra and half-lives 
of zinc isotopes produced through photonuclear reactions. The photonuclear reactions were 
induced by a bremsstrahlung photon beam generated by a linac. The particularly interesting 
point, in which we differed from previous such experiments, is that we have used a clinical 
linac. Our aim was to demonstrate the potential usefulness and the power of such machines, 
building on the pioneering work of Ref. 13j. Our work was motivated by the availability of 


clinical linacs. Their presence is on the increase and after decommissioning from their use 
in medical treatments many are readily available to be used in research. In addition even 
linac still in medical use can be utilized for nuclear physics research, with some planing and 
good coordination with the medical institution operating the linac. Although this concept 
and the idea have been demonstrated in Ref. il3|, it has not created a trend of clinical-linac 
use in science. With this paper, we wished to bring this idea back into focus. Our aim was 
to show that a clinical linac can be a competitive tool in modern nuclear physics, and that 
it can be especially useful for laboratories such as ours, in the developing world, as an easy 
and accessible way of performing experiments. 

We used a beam of 18 MeV endpoint energy, well above the proton and neutron sep¬ 
aration energies of all zinc isotopes, activating all of them. However, our experiment was 
not designed to observe prompt gammas since it was offline. The elements, of which transi¬ 
tions from their nuclear levels were observed, are the decay products of radioactive elements 
created by the photo-activation of zinc. These transitions, as well as the half-lives of the 
parent nuclei, were the intended goal and in this respect the experiment performed was quite 
successful. 

The experiment was separated into several parts, the crucial ones being the sample spec¬ 
trum measurement, analysis and the calibrations. On the calibration side, we have made 
substantial efforts to understand the sources of errors and perform the analysis as accurate 
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and complete as possible. We have performed and combined two calibrations, before and 
after the sample measnrement, in order to acconnt for channel drift during the course of 
the experiment. Such a procedure ensured that our data had a good robustness and was 
reliable. We have paid close attention to the htted polynomial, choosing cubic as the best. 
For the sample spectrum, we paid close attention to any background contribution especially 
when considering the area under the peak, i. e. the counts. The second part of the analysis, 
the study of parent nuclei half-lives, was devoted to this. 


As for the results themselves, we have demonstrated quite conhdently that experiments 
such as ours improve the accuracy of transition energies. Indeed, we have obtained results, 
which are in good agreement with the literature values, but notably often with reduced 
uncertainties. Our main contribution can be seen in the case of ®^Zn, whose levels we have 
consistently determined to an accuracy level that is same or better than the one found in 
the literature. Granted the measurement to which we compared our data came from an 

n 

experiment performed nearly 40 years ago Ref. but this is exactly where we expect to 
improve on data from the literature. The results obtained and the setup of the experiment 
make the study presented in this paper quite interesting and a valuable contribution to the 
nuclear data set. In fact, we hrmly believe that such clinical-linac-based studies on proton- 
rich nuclei can offer improvements of data in many intermediate-mass nuclei. In accordance 
with this idea we are currently investigating the application of the our experimental setup 
to other nuclei. Currently we have already performed several follow up measurements on 
nuclei such as Cl, Br, Sc, Ga, Sb and Pr and are in the process of analyzing the data. The 
initial results show much promise. 


On the half-life side, our data are not as nearly impressive as for the transition energies. 
However, given the limitation of this study, it is still a good check of consistency. Further¬ 
more, with the improvement of our experimental setup, primarily in target preparation and 
choice of nuclei, we are conhdent that further studies will offer better results on the half- 
lives as well. Obviously, increased experience in this kind of experiments will also contribute 
to improvement of results in the future. In fact our follow up experiment on other nuclei 
indicate a good improvement in the half-life determination. 
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